The low-lying energy-momentum spectrum for the four-Fermi model on a lattice 
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We obtain the low-lying energy-momentum spectrum for the imaginary-time lattice four-Fermi or 
Gross-Neveu model in d + 1 space-time dimensions (d = 1,2,3) and with iV-component fermions. 
Let k > be the hopping parameter, A > the four-fermion coupling and M > denote the fermion 
mass; and take s x s spin matrices, s = 2, 4. We work in the k C 1 regime. Our analysis of the one- 
and the two-particle spectrum is based on spectral representation for suitable two- and four-fermion 
correlations. The one-particle energy-momentum spectrum is obtained rigorously and is manifested 
by sN/2 isolated and identical dispersion curves, and the mass of particles has asymptotic value 
— In k. The existence of two-particle bound states above or below the two-particle band depends 
on whether Gaussian domination does hold or does not, respectively. Two-particle bound states 
emerge from solutions to a lattice Bethe-Salpeter equation, in a ladder approximation. Within 
this approximation, the sN(sN/2 — l)/4 identical bound states have O(k ) binding energies at zero 
system momentum and their masses are all equal, with value ~ —2 In k. Our results can be validated 
to the complete model as the Bethe-Salpeter kernel exhibits good decay properties. 

PACS numbers: 11.10.Kk,11.10.St,02.70.Hm,71.10.Fd 



We consider an imaginary-time lattice functional in- 
tegral formulation of the four-Fermi or Gross-Neveu 
(GN) model with A-flavor component fermion fields 
in d + 1 dimensions, d = 1,2,3, and a global U(A) 
flavor symmetry. 

The interest in such a model is because it is very 
rich and presents a non-trivial behavior. For exam- 
ple, in d = 1, the GN model [lj is renormalizable and 
asymptotically free in the ultraviolet (N > 1); it also 
exhibits dynamical symmetry breaking with mass gen- 
eration and dimensional transmutation, and asymp- 
totic completeness has been verified up to the two- 
particle threshold [2|. Also,the GN model is simple 
enough to be, in principle, accessible via analytical ap- 
proaches It has been employed to test properties 
of more complex field theories. It is a purely fermionic 
model and, since fermions are bounded operators, it 
is in a special position to apply (non-)perturbative 
techniques to obtain rigorous control of some correla- 
tions properties, going beyond the level of simply few 
perturbation orders [J]. Moreover, the GN model is 
important for the description of interesting physical 
phenomena as superconductivity (see e.g. Ref. [5|). 

Let k 3* be the hopping parameter, M > the 
fermion bare mass, A > the four fermion coupling 
and take s x s spin matrices. We work in the s<l 
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regime and analyze the one- and two-particle sectors 
of the energy-momentum spectrum. We show that 
there are sN/2 particles manifested by isolated dis- 
persion curves w{p) = w(p, k, M, A), up to near the 
three-particle threshold (upper mass gap property). 
From the U(iV) symmetry these curves are identical; 
the particle masses m = m(«) = w(p = 0) are all 
equal and with asymptotic value — In k. Two-particle 
bound states near the two-particle band are obtained 
from solutions to a lattice Bethe-Salpeter (B-S) equa- 
tion, in a ladder approximation. Their existence above 
or below the two-particle band depends on whether 
Gaussian domination does hold or not, respectively. 
Within the ladder approximation, the sN(sN/2 — l)/A 
bound state masses are all equal and of order — 2 hiK, 
their binding energies are 0(1). The methods applied 
here, were established in Refs. [1, 0, H> H Efll for lattice 
spin and field theoretical models. Also, our method 
has the main ingredients to extend the bound state 
results beyond the ladder approximation (see Refs. 

Bl)- 

Precisely, our system action is given by 

-(A/AT) E x ^a,a(x)] 2 , 

with summation over repeated indices, and where 
ipa,a.(z)) tl>p,b{y) are Grassmann variables, x — (x°,x) 
= (x°,...,x d ) is a lattice site, a; £ Zi/ 2 = 



2 



{±1/2, ±3/2, ...}, x e Z d , a, (3 € {1, 2, .., s} are spin 
indices , a,b £ {1,2,. ..TV} are flavor indices, e = ±1 
and, for /i = 0,...,g?, the e M are unitary vectors of 
Z d+1 . Originally, instead of M, the lattice formulation 
gives us a bare mass M.(fh, k) = (m+2«;) but, without 
loss of generality, given k we choose the bare fermion 
mass m — 0(1) » n such that M.(fh, k) = M > 0. 
Also, Y ±e ' = —7 ± where if are <i + 1 hermi- 
tian traceless ant i- commuting matrices. For d = 2,3, 
we take s — 4 and <y9 M = 7 M , [i = 0,1, ... ,d, are 

4x4 Euclidean Dirac matrices 7 = ( q 2 _° 2 , 

7 J = ( ^ ^j, and cr J , j = 1,2,3 are the her- 

mitian traceless anti-commuting Pauli matrices. In 
these cases, we have of course s = 4. For d = 1, 2, we 
also consider the choice if = a 3 , ipj =1 ----- d = being 
Pauli matrices, and s = 2. From now on, we concen- 
trate our analysis here only on the d = 3, s = 4 case. 
The other cases can be reproduced from this one. 

We start in a finite-volume A C = Z 1/2 <g> Z d . 
The normalized expectation of a function F(ip,ip) is 

Z A J A 

where dip dtp = l\ x a> a dip a . a {x) dip a ^ a (x) and Z A = 

Ja e ~ SI '^'^ S> d-i{; dip. Using a polymer expansion [ll|, we 
show that correlations are given by convergent expan- 
sions if (k/M) <C 1, and that the thermodynamic limit 
A — > Zq +1 exists. The limiting correlation functions 
(cf 's) are denoted by ( • } and are translation invariant 
and analytic functions jointly in k and A. As a conse- 
quence, the truncated cf 's have exponential tree decay 
with rate at least —(1 — e) ln[0(l)/i], where e — > as 
k — ► 0. Below, we work in the thermodynamic limit. 

The choice of the shifted temporal lattice is to ac- 
commodate two-sided equal-time limits of Fermi cor- 
relations in the continuum limit. For n ^ the action 
of Eq. preserves positivity and, therefore, we can 
use a Feynman-Kac (F-K) formula (see [111]) to con- 
struct an imaginary-time lattice quantum field theory 
associated with this model. This construction pro- 
vides a positive metric Hilbert space H, with inner 
product (F,G)n, vacuum vector f2 and field opera- 
tors ip^{x), ^(x). A straightforward consequence of 
this construction is that the squared transfer matrix 
Tq = e~ 2H is, in the infinite-volume limit, a well de- 
fined self-adjoint contraction operator in the Hilbert 
space of physical states TL (assuming that To does not 
have a eigenvalue). Thus, we can construct com- 
muting self-adjoint energy-momentum (E-M) opera- 
tors H ^ and P — (P 1 , P d ) (with spectral points 
p G T = (— tt, n] d ), the generators of time and uni- 
tary space translations T — e lP - x , respectively, and 
analyze the joint spectrum of Tq and P. 



For F and G with fixed time, and for the uppercase 
T denoting an operator on H, taking a: > 0, the F-K 
formula is {G,ffffffffF) H = ([Tq° T s F]QG) , 
where T x = Tf ...TJ and is an anti- linear opera- 
tor which involves time reflection. The key idea of our 
method to detect particles and bound states in the E- 
M spectrum is to use the F-K formula to obtain spec- 
tral representations for conveniently truncated two- 
and four-point cf's. It is these spectral representa- 
tions that allow us to relate the complex momentum 
singularities of the Fourier transform of these cf's to 
the E-M spectrum (see @). To obtain the one-particle 
spectrum, we analyze the two-point cf (x° — y 5^ 0) 

G af}{ X ^y) = Xy«>xO (lpaa{x)lpf3b(y)) ^ ,^ 
-Xy°<x° (^a a (x)lp/3b(y))* , 

where x( x ) is the Heaviside function. Using F-K and 
the spectral representations of To and P, we obtain 
the spectral representation for the two-point cf 

G%{x,y)=f f xf-^e^-y^dp^X), 

(3) 

where dp a ^{X) = d(ipp a ,£{X)ip a b)n with E(X) = 
E o(Xo)U.j =1 F j (X j ), and where T (A ) (Fj(Xj)) are 
the spectral families of e~ H (-P- 7 ). 

To obtain the two-particle spectrum, we analyze the 
partially truncated four-point cf 

D^ 3 (x 1 ,X 2 ,yi,y2) = Xy°^x° 1 ($<xia,i(xi)$a 2 a, 2 (x2) 
X IppbiiVlWfcbzfa)) -Xj,°>x° (^cuaiOl) 
X 1p a2 a 2 {x 2 )lpf3 1 b 1 {yi)lpf3 2 bAy2))* , 

(4) 

where ab stands for a\a 2 b\b2, ot(3 for a,\u 2 pi\Pi 2 , and 
we work in the equal-time representation x\ = x% and 
Vi = 2/2- Using the lattice relative coordinates £ = 
x 2 - xi, ff = 2/2 - yi,r = yi - x 2 (see Ref. 0,0), we 
obtain, as for G, the spectral representation (r° 7^ 0) 

D$&n,T)= f I X^-'e^+^dv^X), 

where dv%{X) = d{^%{- V ),£{X)r a \Z{0)-H and 
<P^(x) = ipl{Q)ip b p{x)Q.. Without spectral represen- 
tations, the analysis of the decay rates of cf's does 
not lead to particles and bound states, specially when 
multiplicities are present with narrow separations. 

The cf's given in Eqs. @ and ^ are matrix valued 
operators on convenient spaces, and we can use sym- 
metries at the level of cf's to simplify them. We follow 
Rcf. 10] where symmetries were analyzed. A symme- 
try operation Y on the Grassmann algebra Q of the 
fields ip and ip is an operation of Q onto itself such that 
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of cf's: y(x°,x) — (— x°,x*) and A = i 



^lpaa(x) = Aapljlpbiyx), Y^ aa (x) = 1ppb(yx)Bi3 ai 

where A and B are 4x4 complex matrices and y is a 
bijective linear map of Zq +1 onto itself. The symbol ~ 
stands for the introduction or the removal of a bar. A 
symmetry of the model is a symmetry operation which 
leaves invariant the action (fT]) and the cf's (eventu- 
ally apart of a complex conjugation operation). Usual 
symmetries as discrete spatial rotations, time reversal 
T , charge conjugation C and parity V can be imple- 
mented as unitary (or anti-unitary for time reversal) 
operators in 7i are symmetries of the model. We also 
found a time reflection symmetry T, only at the level 

' -I a 
h 

B = A -1 . By applying C, the particle and the anti- 
particle cf's are the same, and we only consider lower 
spin indexes, i.e. a, ft = 3, 4 in the cf's of Eqs. © 
and (gD. Next, using CVT, we see that G of Eq. © 
is a multiple of the identity. So, we drop the spin and 
flavor indices (due to the U(iV) symmetry!) from G. 

Concerning the one-particle spectrum, we follow the 
analysis of [6|, for spin models. Using the hyperplane 
decoupling method ([1(1 Q3)> we see that G(x,y) is 
bounded by \G{x,y)\ < 0(1)k^-^. Thus, taking 
G = Gd + G n as a matrix operator on £ 2 (Z,q +1 ), Gd = 
G(x, x) ^ 0, we can define the inverse convolution 
matrix V = G -1 through a convergent Neumann se- 
ries. Using the hyperplane decoupling method, we can 
also show that the improved global bound \T(x,y) = 
T(x - y)\ < 0(1) k 3 ^"^! holds for \x - y\ > 1. Also, 
by a generalized Paley- Wiener theorem, the Fourier 
transform T(p) — J2 X F(x)e~ tp - X is analytic up to close 
to the three-particle threshold « — 31nK. Hence, from 
Eq. ((3]), the one-particle spectrum is determined by 
the zeroes of T(p). By the analytic implicit function 
theorem and a Rouche's theorem argument, there is 
a zero of multiplicity sN/2 (due to spin/flavor) corre- 
sponding to isolated particle dispersion curves 

u)(p) = — hig 2 K — g 2 K cosp J '+ 0(k 2 ) , (5) 
3=1,. ...d 

where g2 = ( - Qa (a;)V'aa(a;)}U=o- The particle masses 
are m = m{n) = w(p = 0) = — lng 2 K — dg 2 n + 
0{k 2 ). Above the one-particle spectrum there are 
two-particle,... finite width free-particle bands which 
eventually overlap for a large particle number. 

To search for two-particle bound states, we use a 
lattice B-S equation for D. In operator form, it reads 

D = D + DKD . (6) 

In terms of its operator kernel, Dq is given by 

D* b a0 {x,y) = G a a \%(x 1 ,y 2 )GZ%(x 2 ,y 1 ) 



G%%(x uyi )G%%(x 2 ,y 2 



a 2 b 2 



and is the Gaussianly evaluated four-point function 
with the corrected propagator. Formally, from Eq. 
©, K = Dq 1 - D- 1 . It is easy to show that D, D 
are bounded operators in £ 2 (A), the antisymmetric 
subspace of £ 2 {A), A = {(x,y) 6 1 d Q +1 \x° = J/T^- 
Next, by a Neumann series argument, we show that 
Dq 1 , D~ x and hence K~ x exist as bounded operators 
in £ 2 (A). This uses decay properties of Dq. 

Computing K to the leading n (ladder) approxima- 
tion K = L, and using this to solve Eq. © and to 
find the singularities of D, we obtain the bound states. 

We now show that L is local and O(k ). For this, we 
write D = D° + D T . Thus, D T is the connected (trun- 
cated) four-point cf, and vanishes at n = for non- 
coincident points (see Ref. @), as seen by writing D T 
in terms of source derivatives of the logarithm of the 
generating function. Using a Taylor expansion at n = 
0, we find D$,(x,y) = -2g A P xy {\ - Sf p )l A -2g 2 (l - 
P xy )l A + O(k) and D^ ap (x,y)) = ~2g 2 l A +0(k), 

where P X y — 8x 1 x 2 8y i y 2 8y- LX2 t = 8aia 2 8/3 1 fi 2 8fi ia2 

8 ai a 2 Sbib 2 8b! a 2J is the identity on 4%(A) and 
9i = {ipaa(x)4!iSa{x)ip/3a(x)ip a a(x))\ K =a. With all this, 
L%(X,V) = -mtr 1 ~ ^9a)- X ]P xv [1% - S%] = 
-XP xy [lf - S^ b ], where I is the AN 2 x AN 2 identity. 

In the ladder approximation, the B-S equation reads 
Df (x,y) = D$ a p{x,y) +Y, z D%-{x,z) Lf e (z,w) 
D^ b e/3 (w,y). Changing to the coordinates £, r\ and 
r, taking the Fourier transform in r only, and set- 
ting /(£,??, = E T /(^,T)e^ 0T °, at zero spatial 
momentum, we find Df (^,f),k°) = P> ab afi (t f), k°) 
-D^ aC (C0,k a ) -S%\ Df (6,if,k ). Therefore, 

we obtain (suppressing k°) ff) = D% b a/3 (£,ff) 

+D a o % c €0) {[I + m - S}D (d,0)}- x }£ Dfcp&ff). 
From the analysis of G, D^ b a p(^,f],k ) is analytic in 
[0,2m(/t)) U (2m(«) + W,3m), where W = 2[w(tt) - 
w(0)] — 8ng 2 d + 0(k 2 ) is the width of the two- 
particle band. Hence, there is no bound state for 
H = 0. Bound states occur only when det{I + H[I — S] 
Dq(0, 0, fc)} — 0, and then are associated with a zero 
eigenvalue of I + K [I - S] D o (0, 0, fc°). By a straight- 
forward calculation, D (0,0,t) = -2G 2 (0,r)(I - 
S), then we just have to obtain the eigenvalues of 
I - 2Ki?(fc°)(I - S), where R(k°) = G 2 (0,fc°) = 
/ (4>ia{0)ipia(T)) 2 e- lk ° T ° dT dr. There are 4A^ 2 eigen- 
values. The eigenvalue Ai = 1 is (2N 2 + A r )-fold 
degenerated, has eigenvectors in the symmetric sub- 
space of R 4JV and does not satisfy the bound state 
condition. The eigenvalue A2 = 1 — 2Hi?(fc°) is 
(2A^ 2 — AT)-fold degenerate and has eigenvectors in 
the antisymmetric subspace of R 4Ar2 13]. These 
states may satisfy the bound-state condition. For 
them, the condition is given by R(k°) = (2H) _1 . 
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We next use the analysis of G to write Eq. ([3]) 
as (a; ^ 0) G(0,x) = JT/ T<I e l ^ s ~ E \ x ^dap(E)dp. 
Also, for in Eq. ([Sf, we use the decomposi- 

tion dap(E) = Z(p, k)S(E - u){p))dE + dap(E), with 

z (p,k) = f|(p° = ix,^)| x=w (^), and where f(p) = 
[G(p)] _1 . Both, dap(E) and dap(E), are positive mea- 
sures. dap(E) has support in (m, oo), m w — 3hiK 
is a lower bound for the onset of three-particle spec- 
trum. The decomposition of dap(E) dp gives a separa- 
tion into one-particle and the three- or more particle 
contributions (given by daff(E)). Thus, Z(p, k) > 0. 

Proceeding as in Ref. [g, Q , and holding only the 
product of one-particle contributions in G 2 (associ- 
ated with Z(p, k)S(E — uj(p))dE), we obtain 



i?(fc°) = 2(27r) 



d+l 



sinh(2o> (p))[Z(p, K )} 2 
■d cosh(2w(p)) — cos/c° 



dp+0(n) 



where Z(p, k) = (2Tr)- d g 2 + 0(K). So, for 0<3(/c°)< 
2m (i.e., below the two-particle threshold), R(k°) > 0. 
Hence, there are no bound states for H < 0. Similarly, 
for 2m + W< Q(fc°) < 3m, R(k°) < and again no 
bound state exists for K > 0. The only k° singularities 
of £>(£, if, k°), for 3f(fc°) £ (0, 2m) U (2m + W, 3m), are 
solutions of 2K£ T G 2 (0,T)e 2mr V Ar ° = 1, and we 
have set ik° = —2m + A (recall m is the mass). 

We next give an intuitive argument for the bound 
state formula based on the behavior of (^i a (O) , 0ia(''"))- 
A rigorous argument follows from using the convo- 
lution form in momentum space of the above con- 
dition and Eq. (J3j) . Expanding (i'i a (0)ipia(T)) to 
the leading order in k, we obtain (4'ia(0)' l Pia( T )) ~ 
K \ T l+l 1 "!^ 1 " ' + ' r ' +1 . This follows from expanding the 
numerator of the two-point cf. For the fermion inte- 
gration to be nonzero, there must be a chain of over- 
lapping bonds connecting the two points. Holding the 
chain of minimal length, the above result follows. A 
polymer expansion is used to rigourously control all 



contributions. Using the value of m, we obtain the 
bound state condition, to the leading order in k, as 
2g 2 X(l-e- A ) = 1. 

Let A A = A (A R = -A - W) denote the binding 
energy at zero system momentum for the attractive 
(repulsive) case. For the attractive case, we find A A = 
-ln[l - (2^g 2 y 1 } : where N = [g 4 - g 2 2 ]/(2 3 2 34 ) > 
(Gaussian Subjugation, as the four-point function 
dominates the square of the two-point function), giv- 
ing a bound state below the two-particle band. For 
the repulsive case, A R = ln[l — (2K52) 1 ] — W, where 
now K < (Gaussian Domination), and the bound 
state is above the band. 

In closing, we point out that the low-lying E-M 
spectra of other GN models in d = 1 have been ob- 
tained (at least approximately) by other methods in 
the continuum massless case. For instance, the Bethe- 
Yang ansatz was employed in the chiral invariant ver- 
sion of the GN hamiltonian, as well as semiclassical 
approaches as in Ref. Q. Our results agree qual- 
itatively with those. However, besides the fact they 
are obtained in a precise mathematical setting and in- 
clude the upper gap property (absent in hamiltonian 
formulations), and that our method can be easily ex- 
ported to other models, our results can be extended 
beyond the ladder approximation using the methods 
e.g. of Refs. [8, 9]. To do this, a fast decay of the B-S 
kernel is crucial. Using the hyperplane method, we 
get \K$,(£,t},t)\ < O(l) K 4 l- l+l(ll«+?+ 2? ll+ll«ll + ll^ll)|. 
Whether or not the degeneracy of the bound states is 
preserved in the complete model is a good question. 
Finally, we remark that the spectral pattern obtained 
here for the U(iV) GN model is similar to the one 
found for the 0(2A^) spin model in Ref. 0], except 
that the later has also a different bound state of mul- 
tiplicity 2JV. 
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